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For a Gaussian spectrum of primordial density fluctuations, ultracompact minihalos (UCMHs) 
of dark matter are expected to be produced in much greater abundance than, e.g., primordial 
black holes. Forming shortly after matter-radiation equality, these objects would develop very 
dense and spiky dark matter profiles. In the standard scenario where dark matter consists of 
thermally-produced, weakly-interacting massive particles, UCMHs could thus appear as highly lu- 
minous gamma-ray sources, or leave an imprint in the cosmic microwave background by changing 
the reionisation history of the Universe. We derive corresponding limits on the cosmic abundance of 
UCMHs at different epochs, and translate them into constraints on the primordial power spectrum. 
We find the resulting constraints to be quite severe, especially at length scales much smaller than 
what can be directly probed by the cosmic microwave background or large-scale structure observa- 
tions. We use our results to provide an updated compilation of the best available constraints on 
the power of density fluctuations on all scales, ranging from the present-day horizon to scales more 
than 20 orders of magnitude smaller. 
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I. INTRODUCTION 

Structure formation in the Universe is thought to have 
proceeded in a bottom-up process, starting with the for- 
mation of smaller bound objects from the gravitational 
collapse of small-amplitude, random initial density per- 
turbations; in a second step, these smaller objects would 
then merge to larger ones. In the simplest scenario [IHS] 
(see also Refs. [3 |S] for an introduction to the standard 
model of cosmology), a nearly scale-invariant (Harrison- 
Zel'dovich) [H [10] Gaussian spectrum of perturbations 
is produced during inflation, and the bulk of perturba- 
tions are of sufficiently small amplitude at horizon entry 
{6 ~ 10~^) that they do not collapse until well after 
matter-radiation equality. 

None of these assumptions need be fulfilled, however. 
Many inflationary models (and some alternatives) predict 
a departure from scale invariance (see e.g. Refs. pTMlQj ). 
or from Gaussian statistics (see e.g. Refs. [^OH^ and ref- 
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erences therein). In theories with more complicated cos- 
mological histories, events such as phase transitions and 
late-time dynamics of scalar fields such as string mod- 
uli, the inflaton or curvaton [251 - 130] might have lead to 
an injection of additional power on specific scales. Re- 
gardless of when the primordial spectrum was produced, 
the presence of additional power on some scale may have 
led structures of a size corresponding to that scale to 
collapse far earlier than in the canonical scenario. In se- 
vere cases, this collapse may even have occurred before 
matter-radiation equality |31j . 

The most extreme and best-studied example of such 
rapid gravitational collapse is that of primordial black 
holes (PBHs) A PBH is expected to form when 

a perturbation enters the horizon with such a large am- 
plitude {5 > 0.3 — 0.7) that a substantial fraction of the 
horizon volume collapses directly to a black hole [351 [SS] ■ 
Even very large density perturbations with 5 > 1 form 
PBHs, as has only recently been clarified 'Tf\, rather 
than closing up upon themselves and forming separate 
universes, as originally argued by Hawking |33j . PBHs 
have been proposed as possible dark matter (DM) can- 
didates j38.- 41j . but the limits on their abundance are 
tight [111 |43], and the formation process lacks a natu- 
ral mechanism for producing the observed cosmological 
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abundance of DM. Furthermore, due to the enormous 
amphtudes required for their formation, a sizable PBH 
abundance can only be produced in scenarios that devi- 
ate rather strongly from the simple Harrison-Zel'dovich 
spectrum ^40t .44.-46, . 

Less severe departures, at the level of (5 ~ 10"'^ dur- 
ing early radiation domination, can instead lead to the 
formation of so-called ultracompact primordial miniha- 
los (UCMHs) gZHin]- Such a perturbation in the DM 
component will continue to grow until it eventually un- 
dergoes gravitational collapse before or very shortly af- 
ter recombination. The mass of these objects is propor- 
tional to the horizon mass at the time such a density 
perturbation enters the horizon, i.e. it is uniquely deter- 
mined by the wavelength of the perturbation mode. In 
general, smaller-scale perturbations of a given amplitude 
will more easily produce UCMHs, as they enter the hori- 
zon earlier and thus give the over-density more time to 
grow. Unlike in the PBH case, the lack of an event hori- 
zon means that the seed cannot retain radiation (which, 
at this early time, includes baryons), so baryonic mat- 
ter and photons free-stream immediately back out of the 
over-density. The seed thus consists entirely of cold dark 
matter (CDM), but it will later accrete both dark and 
(after recombination) baryonic matter. 

The amplitudes of temperature inhomogeneities in the 
cosmic microwave background (CMB) [SHI [ST] allow us to 
infer the power of the primordial density fluctuations on 
scales that range from 0(10* Mpc), almost the horizon 
size today, down to 0(10 Mpc), roughly an order of mag- 
nitude smaller than the size of the sound horizon at the 
time of recombination. The observed large-scale distribu- 
tion of matter in the Universe also carries clues as to the 
amplitude of linear cosmological density perturbations, 
on scales from about 10"^ Mpc down to scales somewhat 
smaller than what is currently accessible by CMB obser- 
vations [52] ■ Weak gravitational lensing gives access to 
scales down to 1 Mpc [S3j while measurements of the in- 
tergalactic hydrogen clumping, traced by the Lyman-a 
forest, constrain the density perturbations on scales as 
small as O(10~^)Mpc [S3J[S3]. Observational evidence 
for the presence of cosmological perturbations at smaller 
scales is essentially absent; the best current limits come 
from the non-detection of PBHs ^43] HI]. Given their 
sensitivity to the amplitude of small-scale perturbations, 
UCMHs are potentially much stronger probes of the pri- 
mordial spectrum at large k than PBHs [13 ES] • 

The most popular and compelling model for DM is 
that it consists of weakly-interacting massive particles 
( WIMPs) , thermally produced in the early Universe [SSj- 
In this case, the compactness of UCMHs makes them 
prime targets for indirect detection of self-annihilating 
DM; the non-observation of corresponding gamma rays, 
in particular, places constraints on the UCMH mass frac- 
tion. The most sensitive searches for gamma rays from 
DM annihilation have been with Air Cerenkov Telescopes 
(ACTs) [BUHSg and the Large Area Telescope (LAT) [63- 
[7T| aboard the Fermi satellite |72j . Another promising 



technique is to consider the contribution of DM anni- 
hilation within UCMHs to the reionisation of the inter- 
galactic medium at redshifts z ^ 6 |73j . By comparing 
the predicted Thomson-scattering optical depth from the 
present day to the surface of last scattering, Tq, with the 
value observed in the CMB, UCMH fractions leading to 
very early reionisation (and therefore large r^) can be 
excluded. 

In this article, we provide detailed and updated limits 
on the primordial spectrum at small scales. In particu- 
lar, we include the non-observation of UCMHs to date 
by gamma-ray searches for DM with the Fermi-h AT , 
as well as potential impacts of UCMHs upon reionisa- 
tion. We expand on previous works in this direction, 
most notably Ref. [55) . by deriving from first princi- 
ples the minimum density contrast (S™'" required to form 
a UCMH, discussing in detail the transfer and window 
functions needed to correctly normalise the primordial 
power spectrum to the one observed today, and by using 
an improved treatment of the statistics of Fermi non- 
detection of individual sources. For comparison, we also 
compile and incorporate other constraints on the primor- 
dial power spectrum into our final limits, extending from 
the horizon size today down to scales more than 20 orders 
of magnitude smaller. 

We begin in Sec. [TTj by describing the formation and 
structure of UCMHs. In Sec. |III[ we then briefiy recap the 
calculation of their cosmological abundance and provide 
limits on their present number density. We continue in 
Sec. |IV| by translating our limits to new constraints on the 
primordial power spectrum from UCMHs and discussing 
existing ones, before concluding with Sec . [Vj We give a 
detailed calculation of 5™™ in Appendix lAl and of the 
mass variance of perturbations in Appenmx [B] taking 
into account the correct window and transfer functions. 



II. FORMATION AND STRUCTURE OF 
UCMHS 

As noted in Ref. (48j, density contrasts as small as 
6 ~ 10^'^ during radiation domination already suffice to 
create over-dense regions that would later collapse into 
UCMHs (see Appendix [A| for more details). The initial 
mass in CDM particles contained in such a region, i.e. the 
mass at the time the corresponding fluctuation of co- 
moving size R enters into the horizon, is given by 



J aH=l/R 



Vo.ii2y VMpcy °' ^ ' 

where p^^ refers to the actual DM density at the time of 
horizon entry, and fi^ is the fraction of the critical den- 
sity of the Universe in DM today. Note that we define R 
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here as the co-moving radius of the collapsing region; the 
co-moving diameter of the region is 2i?, which is half the 
co-moving wavelength of the corresponding physical den- 
sity fluctuation. No factors of gcs, the effective number 
of relativistic degrees of freedom at various times in the 
early history of the Universe, enter in the above expres- 
sion, because only the CDM component collapses and 
contributes to the initial UCMH mass. 

During radiation domination, the UCMH-forming 
'seed' consists only of CDM, and its mass stays essentially 
constant. Around matter-radiation equality, it then be- 
gins to grow by infall of both dark and (after decoupling) 
baryonic matter as 

1 



UCMH 



(^) 



z + l 



(2) 



We conservatively assume that this growth only continues 
until standard structure formation has progressed suffi- 
ciently far as to allow star formation, at z ^ 10. After 
this time, dynamical friction between DM halos and hier- 
archical structure formation presumably make further ac- 
cretion from the smooth cosmological background of DM 
inefficient (note that the growth rate presented in Eq. ^ 
assumes accretion in a homogeneous and unbound DM 
background). Taking fi^^/i^ = 0.112 [5T], the mass of a 
UCMH today is then related to the size of the original, 
slightly over-dense region by 



UCMH 



UCMH 



(z < 10) w 4 X 10 
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Similarly, the current-day UCMH mass is related to the 
horizon mass Affj {R) at the time when the fluctuation of 
co-moving size R enters the horizon, by^ 
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where ^cff (Soff) evaluated at the time of horizon entry 
(equality) . 

After kinetic decoupling |74] , the velocity of CDM par- 
ticles decreases as w oc (1 -I- z); the DM within the over- 
dense region therefore initially has an extremely low ve- 
locity dispersion as long as the density fluctuations are so 
small that they do not induce sizable gravitational poten- 
tials. Even shortly after the onset of matter-domination, 
when the fluctuations that we are interested in here start 
to become non-linear, the velocity dispersion increases 
only mildly as |48 t [75 ] 
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1 Eq. |4| is rigorously correct only when R <C (acqHcq) ^ ~ 1.0 X 
10^ Mpc, i.e. Mh < M'j^ ~ 3.5 x lO^Mg; close to Mh = M'^^ 
it receives a further correction of 2~^/*. We give this expression 
here only for explanatory purposes, and use Eq. which is 
exact, for all the calculations that follow. 



with av,o = 0.14 ms""'^. UCMHs thus form by almost 
pure radial infall [48] . which leads to the growth rate 
presented in Eq. ([2|, and a DM radial profile 



/ N 3/xMucmh(z) 
p^[r,z) ~ 3—; 

167ri?uCMH(^)*''' 
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where = Q^/flm denotes the fraction of matter that 
is CDM. This extremely steep profile, a direct conse- 
quence of spherically-symmetric collapse, was first de- 
rived analytically [7S1 [77] and later confirmed by explicit 
A^-body simulations [75] [75] (see in particular Fig. 6 of 
Ref. [78j). The factor enters because even though the 
initial UCMH seed consists only of DM, the matter it ac- 
cretes following decoupling includes both dark and bary- 
onic matter, and accounts for the majority of the mass of 
the UCMH at the present time. R\jcmh{z) refers to the 
effective radius of the UCMH at redshift z, beyond which 
the density contrast associated with the (fully collapsed) 
UCMH is S < 2. For coUisionless DM, this turns out to 
be |77j 



^ucmh(z) ~ 0.339i?ta(z), 



(7) 



where the turnaround radius i?ta(2) is the radius within 
which matter contained in a collapsing perturbation sep- 
arates from the Hubble flow. Rta,{z) has been obtained by 
Ricotti [801 from fits to numerical simulations of matter 
accretion at early times, and then converted by Ricotti, 
Ostriker & Mack [7S] to -Rucmh(2) using Eq. (|7|, giving 



^ucmh(^) 
pc 



With our assumed cutoff in accretion at the beginning 
of star formation at z ^ 10, the current profiles and radii 
are obtained by choosing z '--^ 10 in Eqs. (6j8|. Choosing 
e.g. z ~ 30 instead has a minimal impact on results [49) . 

The steep density profile presented in Eq. ^ is valid 
so long as the infalling dark matter follows an approxi- 
mately radial path. This approximation breaks down in 
the innermost parts of the minihalo, where the average 
tangential velocity Vrot of infalling material exceeds the 
local Keplerian velocity uxcp- This occurs only in the 
inner region because Uiot rises more steeply with decreas- 
ing radius than WKop- The radial infall approximation is 
hence violated at steadily larger radii as time goes on, 
as the velocity dispersion of infalling matter increases 
with time, so matter accreted at later times begins its 
infall with larger Wrot- This tends to suppress the con- 
tribution of newly-accreted matter to the inner parts of 
the halo. Accretion following UCMH formation therefore 
contributes preferentially to the outer parts of the halo 
[151 [77] . but leaves the steep inner profile (established 
by radial infall during the earliest stages of formation) 
essentially intact. 

Even for the earliest-accreted material however, radial 
infall cannot be valid all the way to r = 0. This means 
that the DM profile in Eq. ^ can only be expected to 
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be valid down to some cutoff radius r^in- To estimate 
^■111111, we calculate the Keplerian and average tangential 
velocities of the earliest-accreted material in the UCMH 
as a function of the radial distance r, and then simply 
solve for the value of r at which these velocities are equal. 
With M{r) the total mass contained within radius r, 
from Eq. (|6| we have 



WKep(r) 



GM(r) 
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which we see, from Eqs. ^ and ([s]), is independent of 
redshift. By angular momentum conservation, the mean 
tangential velocity of the infalling gas is 

Vrot{r, z) = . (10) 



Setting UKcp and Wrot equal and using Eqs. (5|8), we find 
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This gives 
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where i?uCMH ^^e present-day UCMH radius. 



(12) 
The 



appropriate redshift to choose here is that at which the 
UCMH in question collapses, z — Zc- This sets the ini- 
tial value of the tangential velocity of material collapsing 
from a height i?ucMH(2c) equal to the typical DM veloc- 
ity CTv(zc) at the time of collapse. As R\]cym{zc) is the 
maximum height from which matter falls into the UCMH 
during the initial collapse, this formalism leads to a con- 
servative estimate of rmin. For a maximally conservative 
estimate of r^^^, one then invokes the largest possible 
value of av{zc) at collapse by choosing the smallest al- 
lowed redshift of collapse; in our case, this is Zc = 1000. 

For r < rmiu, clearly, the central cusp must be flattened 
to some extent, as the r^i profile created by radial in- 
fall must soften as that approximation breaks down. The 
most conservative estimate, which we will follow here, is 
to simply truncate the profile within r^m, resulting in a 
cored profile with central density p = Pxifmin)- Softening 
the central density profile also has the impact of modi- 
fying Eq. ([9]) and hence Eqs. (Ill and ( [T2| . In the most 
pessimistic case, where the density profile is simply trun- 
cated inwards of r„iin, the right-hand side of Eqs. (11 
12) should be multiplied by 2 7 = 2.21. In the limit 
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this results in a reduction of the total 
gamma-ray flux by a factor of 3.28, and a correspond- 
ing increase in the inferred maximum allowed number 
of UCMHs (see Sec. Wv\ by about a factor of 6. When 



translated to limits on the primordial power spectrum, 
such a factor becomes negligible. 

The expression in Eq. ^ for the mean DM velocity 
dispersion is based on a power-law fit to numerical cal- 
culations of CTv as a function of redshift, within spheres 
of different co-moving radii 75 . For each redshift and 
co-moving radius tq, the velocity dispersion is obtained 
by performing the integral 



r{z,roy' 



H{zf 



[l^W^^{kr^)]Vs,{k,z) 



dk 



(13) 

where Vs^ (fc, z) is the power spectrum of DM den- 
sity perturbations. This was computed by the authors 
of Ref. [75 assuming a Harrison- Zel'dovich spectrum 
and a standard ACDM cosmology, using earlier code 
from Ref. |77] to track the coupled evolution of the 
dark and baryonic perturbations. ^ Wth is the Fourier- 
transformed top-hat window function, described in more 
detail following Eq. (B3). In order to obtain Eq. ([5]), the 
fit is evaluated at a co- moving radius equal to i?ucMH(^), 
c.f. Eq. (|8|. The fits have a claimed accuracy of 5%, so 
the impact of the uncertainty in Eq. ([5]) upon our final 
power spectrum upper limits is negligible for power spec- 
tra similar to the Harrison-Zel'dovich case. In principle 
however, the DM velocity dispersion could deviate from 
the behaviour of Eq. ([5| if the calculations of [75 were 
repeated with an alternative input primordial spectrum. 
Investigating this effect is beyond the scope of this paper, 
but might make for interesting future study. 

Finally, one other consideration must be taken into ac- 
count when calculating DM densities in the inner region 
of UCMHs: over time, WIMP DM will annihilate away, 
softening the central density cusp. Following Ref. [5T| 
and earlier UCMH work [33 ES]) we estimate p^,max, the 
maximum possible remaining density at time t in a halo 
born at ij, as 



Px{r < ^cut) = P^, 



{nv){t-uy 



(14) 



Here is the mass of the WIMP and {av) is its velocity- 
weighted, thermally-averaged annihilation cross section, 
taken in the zero- velocity limit. For UCMHs seen today, 
t = 13.76 Gyr [5T]; ignoring gravitational contraction of 
the DM due to baryonic collapse (seen to have mini- 
mal impact upon overall gamma-ray fluxes from UCMHs 
gS]), ti = t{zeq) = 59Myr [SI]. We find that r^t > 4rmi„ 
for all the UCMHs we consider in this paper, so their cen- 
tral densities do indeed violate the annihilation bound, 



^ Vg^ is defined in the analogous way to the total matter power 
spectrum Vs, see Eq. (|B3|, but with the replacement of the total 



matter perturbation 6 by the DM perturbation 5^. Similarly, 
it can be obtained in terms of the spectrum of curvature per- 
turbations V-jz simply by replacing in Eq. | |B22| | the radiation 
transfer function Tr, Eq. ||A9[|, with the DM transfer function 



T^, Eq. (AlOl 
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even taking into account departures from radial infall.^ 
To account for this, we simply truncate our density pro- 
files instead at r = r^uty setting the density within this 
radius equal to p^,max- 

Whilst both are important for total annihilation rates, 
neither the annihilation cutoff nor the correction for vio- 
lation of radial infall have any real bearing upon the inte- 
grated mass of UCMHs. Correcting e.g. Eq. ([s]) to retain 
exactly the same integrated UCMH mass after complete 
truncation of the density profile inwards of r^nin or rcut 
would result in changes at less than the percent level. 



III. THE UCMH POPULATION TODAY 

A. Present abundance 

For perturbations following a Gaussian distribution, 
the probability that a region of co-moving size R (at the 
time when this scale enters the horizon, i.e. 1/R = aH) 
will later collapse into a UCMH is given by 



/3(i?) = 



exp 



d(5. 



(15) 

Here, tr^ jj is the CDM mass variance at horizon entry 
and is the density contrast in the CDM component 
only; in Appendix |B] we provide a detailed recipe for 
computing cr^ ^ . The minimal 6^ required for the col- 
lapse eventually to happen is given by 5™"^ (see Appendix 
\K\ for a precise determination of this value as a function 
of R) and an initial density contrast higher than 6'^^^ 
would lead to the formation of a PBH rather than a 
UCMH. PBH formation is expected above 5™'"'' ^1/4 
from semi-analytical arguments |35| or slightly higher 
values, S'^'^^ ~ 0.5, from numerical simulations [36] (re- 
call that for super-horizon adiabatic fluctuations we have 
~ (3/4) <5 during radiation domination). UCMH for- 
mation, on the other hand, only requires much smaller 
density contrasts of the order of (5™™ ~ 10"^ [481. With 
fv.H ~ 10~^, as observed on large scales, we would thus 



always have ct^.h < < <5y '"^ and thus 



exp 



^min 2 
2<H(i?) 



(16) 



This turns out to be a very good approximation to 
Eq. (15) even in all cases that we will be interested in 



The ratio rcut/fmin reduces as the UCMH mass decreases, so in 
UCMHs Hghter than the smallest we consider here (~10~® ^q) 
the two radii are comparable. The crossing point rcut = ^mln 
occurs at M^Qy^-^ ~ 10~ Mq for the parameters we choose 
here; higher annihilation cross sections or lower DM masses will 
cause this point to shift to higher M^^^j^. 



here, where the power on small scales is significantly 
larger than on large scales. 

Note that /3{R) also counts those regions of size R that 
are contained within a larger region R' > R that satisfies 
5 > (5™'", too. On the other hand, it does not take into 
account the possibility that we have S < J""" for the 
smaller region R, but still have S > i5™'" for the larger 
region R' - in which case the original region of course 
would also collapse eventually and end up in a (bigger) 
UCMH. In the following, we will conservatively neglect 
these contributions to the total UCMH abundance; in 
passing, however, we note that in the Press-Schechter 
formalism |83] these effects would (somewhat arbitrarily) 
be accounted for by multiplying the above expressions for 
/? by a factor of 2. 

Taking into account the accretion of mass described by 
Eq. ([2| , the present density of UCMHs with mass equal 
to or greater than M^^^^ is therefore given by 



^^ucmh(-Mucmh) — ^; 



A/0 



UCMH 



P{R). 



(17) 



where R = Ri^ijcMYi) follows from Eq. Note that 
this expression does not take into account the potential 
destruction of UCMHs due to tidal forces and mergers 
during structure and galaxy formation. Similar to the 
case of super-dense clumps that already collapsed during 
radiation domination [3P, however, these effects turn out 
to be completely negligible. This is because UCMHs form 
so early that by the time of structure formation, they 
have collapsed into quite extreme over-densities with re- 
spect to the smooth background. A good indicator of 
survival probability is the size of the core radius [HH [SS] , 
given by Eqs. (12) and (14). A smaller core radius in- 
dicates a higher survival probability. In particular, a 
core-to-outer radius ratio of less than ~10~^ [5S] indi- 
cates a survival probability very close to unity; for all 
the UCMHs we consider in this paper, this ratio is less 
than 10"*^. It is also worth recalling that UCMHs evolve 
as completely isolated objects for some time after they 
have collapsed: the limits that we will place correspond 
to rather rare fluctuations with 5 /a^^n ~ 3—6 (relative to 
a perturbation spectrum where (Jx,h is already enhanced 
by roughly one order of magnitude on the scale of inter- 
est, compared to what is expected from observations at 
large scales). 

We point out that the DM annihilation signal from 
UCMHs is almost exclusively sensitive to the density 
in the innermost region; even if UCMHs were to lose 
part of their outer material due to tidal stripping, this 
would therefore not affect the corresponding limits that 

In fact, 



we derive in Section IV 
clumps 



even for ordinary DM 
formed in the presence of a standard Harrison- 
Zel'dovich spectrum of density fluctuations and thus 
with much smaller densities than UCMHs - a dense in- 
ner core should remain more or less intact and the im- 
pact of clump destruction on indirect detection prospects 
could be much less severe than one naively might ex- 
pect [551 IHS]; the impact on UCMHs should be even 
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less. For the following discussion, we thus assume that 
Eq. ([3]) indeed provides a very good estimate for the 
present UCMH mass, and Eq. (171 accurately represents 
the present UCMH density. 

As another important consequence of the extremely 
high density of UCMHs discussed above, we note that 
the spatial distribution of UCMHs is expected to track 
the bulk DM. This is quite different from ordinary DM 
subhalos, which are subject to tidal disruption and there- 
fore generally much less abundant in the centres than 
outer parts of large halos, relative to the smooth DM 
component (see e.g. Ref. [57] for a detailed discussion 
and references). Similarly, we expect the effects of the 
stellar disk of our Galaxy on the UCMH distribution to 
be negligible (again in contrast to its effect on ordinary 
DM clumps, which can be sizable [55]). 



B. Limits from gamma rays: individual Galactic 
sources 

If UCMHs consist of WIMP DM, they are generically 
expected to be sources of high-energy gamma rays'' ( [49] ; 
see e.g. Ref. [SO] for an overview of gamma-ray yields from 
WIMP annihilation). In this case, there exists a unique 
distance dobs out to which UCMHs of any given mass will 
be observable by gamma-ray telescopes, given a certain 
instrumental sensitivity. Using the all-sky gamma-ray 
survey performed by the Fermi-L AT , the present abun- 
dance of UCMHs in our own Galaxy can be constrained, 
based on the non-observation of unassociated point and 
extended sources with a spectral signature resembling 
DM annihilation. 

Indeed, no unassociated point or extended sources 
showing evidence of DM annihilation have yet been dis- 
covered by Fermi [TOfM] ; this is true both for sources in 
[nmSg and outside [ST] the 1-year LAT catalogue [95]. 
Although Buckley & Hooper [32] place a rough upper 
limit of 20-60 on the number of DM halos observed in 
the 1 year catalogue, many of these can in fact be associ- 
ated with astrophysical sources; Zechlin et al. ^55] found 
12 possibilities, and then identified the most promising 
as probably a blazar. A recent search in the 2-year cata- 
logue 97] found just 9 potential sources. Given that we 
do not expect all these 9 sources to have been detected 
at better than Scr in 1 year of data, and that statistically, 
we expect at least ~80% of unidentified Fermi sources 
to be relatively easily matched with known sources [55] . 
the implied maximum number of UCMHs in the 1-year 
data is of the order of one or two. Whilst yet to provide 
a statistical upper limit on the number of DM halos, the 
Fermi-JjAT Collaboration itself reports having seen ex- 
actly zero [nH ES llOOj . We thus work under the assump- 



Note that this is true even in the somewhat contrived situation 
where WIMPs annihilate, at tree level, only into neutrinos |89l . 



tion that Fermi observed exactly zero UCMHs during its 
first year of operation, to within its instrumental sensi- 
tivity. Whilst a more detailed treatment would actually 
include a full spectral analysis of all unassociated sources 
in the Fermi survey, such a procedure is well beyond the 
scope of this paper (for a full multi-wavelength approach, 
see Ref. [M])- 

The LAT sensitivity to point sources after 1 year 
of observations, based on a spectral integration above 
100 MeV, is 4 X 10"^ photons cm'^ s'^ [101] for a 5cr de- 
tection. Although this sensitivity is based on a power-law 
spectral source with index —2, expected DM annihilation 
spectra are often sufficiently similar to this that the sen- 
sitivity should be broadly similar. Going beyond this ap- 
proximation would also require spectral analysis beyond 
this paper's scope. We note, however, that pronounced 
spectral features at high energies, close to the DM par- 
ticle's mass (in particular from internal bremsstrahlung 
jl02p . can in principle enhance the effective sensitivity 
by up to an order of magnitude |1031 1104] . 

In order to derive limits upon the fraction / of Galactic 
DM contained in UCMHs, let us for simplicity assume 
that all UCMHs have the same mass -^xjcmh ~ 
sumption we will later comment on. We now pick one 
particular UCMH in the Milky Way, residing some dis- 
tance d from Earth. Assuming that UCMHs track the 
bulk DM, the probability that this UCMH can be found 
within a distance dobs of Earth is 



P, 



M. 



d<c!obs,l 



d<da 



M, 



MW 



(18) 



where Mmw is the total mass of DM in the Milky Way, 
and M(i<d„b^ < Mmw is the mass of DM within dobs 
of Earth. This probability is simply the fraction of the 
(dark) Milky Way mass available for the UCMH to turn 
up in by chance. 

The probability of there existing i such UCMHs within 
dobs can then be constructed from the binomial probabil- 
ity of there being a single one, as done in e.g. Ref. |105j 
for intermediate mass black holes. With the total num- 
ber of UCMHs of mass AfycMH in the MW denoted by 
iViviw, we then have 



P, 



{^^^){Pd<d....ni~Pd<d...,r---^. 

(19) 

Because we assume that all UCMHs have the same mass, 
we can write 



TV, 



MW 



/ 



Mmw 

JX-'"UCMH 



(20) 



where we use Eq. (17) to express the local UCMH mass 



fraction / in the Milky Way as 



/ = r!ucMH/f^„. = /3(i?)/x^^^§^. (21) 

In general, the probability that the number of UCMHs 
i present within dobs is equal to or greater than some 
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threshold number j, i.e. the probabihty that j or more 
UCMHs exist within dobs, is simply 1 minus the individ- 
ual probabilities of there being 0, 1, 2, . . . , j — 2 or j — 1 
of them, 



For a spherically-symmetric DM halo appearing as a 
point source at distance d, the observed gamma-ray flux 
integrated above some threshold energy E^th is 



P. 



d<dabs,i>j 



1 - y^-Pd<rfob3,fc- 



(22) 



k=Q 



So in particular, the probability of there being one or 
more UCMHs within dobs is 



d<da 



,i>l — 1 — PdKdabsfi 

Md<d^ 



1 



1 



MW 



/aJ"mw 



(23) 



This equation provides an estimate of the probability of 
there being any UCMHs at all within an observable dis- 
tance, as a function of the UCMH mass, their maximum 
observable distance, and the UCMH mass fraction. If we 
know the sensitivity for a source detection at confidence 
level (CL) x, then there is a probability x that a UCMH 
residing at exactly dobs, with a mean flux equal to the sen- 
sitivity, would be seen by Fermi. The total probability 
that we would observe one or more UCMHs with Fermi 
is then Ptot = xPd.cd„t,^,i>i- To derive an upper limit 
on / from a non-observation of UCMHs by Fermi, such 
that / < /max at some CL y, we therefore require that 

^tot 

> y for / > /max- That is, y = xPd<d„^^, ' 
which gives 



l/=/max' 



fxK 



UCMH 



log(l-f) 



Ml 



MW 



log 1 - 



Ma 



Mi, 



(24) 



For close objects, M, 



d<dab 



< Ml 



MW 



and we have 



fxKc 



'X^"UCMH ( 1 _ ^ 



(25) 



which is independent of Mmw, as expected. This expres- 
sion has a rather intuitive meaning: it gives the mass 
fraction a single UCMH would contribute to the total 
DM mass within dobs, times a statistical weighting func- 
tion, which accounts for the fact that a non-observation 
of UCMHs does not guarantee that there indeed are none 
in the observable volume, nor does it guarantee that such 
a volume would always contain zero such minihalos, even 
were we certain that this particular volume does not. 

In terms of the gamma-ray flux T{d) from a single 
UCMH at some distance d, this becomes 



fn 



V . 
mill 



UCMH 



log 1- 



(26) 



where py, is the local DM density and the point-source 
sensitivity is defined as J^min = -^(dobs)- Note that this 
expression is in fact independent of d, as T{d) oc d~^. 



nd) = E / 

I. J Eth 



diVfc^^ (akv) 



^p^{r)Ar 



Here {cTkv) is the cross section of the kt\i annihilation 
channel and dNk/dE is its differential photon yield. In- 
verting this expression gives dobs as a function of J^min- 

In order to calculate Md<,da-b!i^ integrate the smooth 
DM profile of the Milky Way over a sphere of radius 
dobs around the Sun, using the NFW profile of Ref. [106] 
(c = 18, M200 = 9.4 X 10" Mq, r, = 17.0 kpc). In prin- 
ciple, local DM substructure in the solar vicinity could 
make this a biased estimator of Md<:daba ^'^^ small dobs- 
Current understanding of the granularity of DM halos un- 
fortunately does not allow the magnitude of this effect to 
be accurately assessed. Recent N-body simulations |107j 
indicate a ^10% variation in density over a sphere of 
500 pc around the Sun, and more over somewhat smaller 
spheres. These scales are at the limit of the simulation's 
resolution however, so going to smaller spheres also re- 
quires higher resolution, which in turn reduces the ob- 
served variation (see the discussion in Sec. 3.2 and the 
Appendices of Ref. |107) ). It is not obvious which effect 
would win out at the smallest scales we consider here, 
which are another 6 orders of magnitude smaller. 



We have constructed Eqs. ( 18 26 1 assuming that all 



UCMHs in the Milky Way have the same mass. Because 
we have observed exactly zero UCMHs of any mass how- 
ever, the resulting limits on / in fact immediately gener- 
alise to arbitrary UCMH mass spectra. This is because 
we know exactly how many UCMHs of each mass have 
been observed (none), so the limits for different masses 
can all be applied independently. If we were instead us- 
ing some j observed UCMH candidate sources to draw 
an upper limit on / (i.e. rather than using them to claim 
a UCMH detection) , we would have a form of Eq. ( 23 ) 



where i > j for some j > 2. In this case, we would have 
to specify the form of the UCMH mass spectrum in order 
to know what fraction of observed sources to attribute to 
each mass band considered in the analysis. Even in this 
case though, the most conservative limits would be ob- 
tained by assuming, for each band individually, that all 
the observed sources were UCMHs with masses in that 
band. 



Using Eq. ( 24 ) and the integrated LAT sensitivity dis- 



cussed above, we have determined the 1-year, 95% CL 
Fermi-LAT upper limits on /, the fraction of Galactic 
DM contained in UCMHs of each mass. This limit is 
shown in Fig. [T] as a solid red line. Here we have calcu- 
lated integrated gamma-ray fiuxes above photon energies 
of 100 MeV, using a suitably extended version of Dark- 
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FIG. 1: The maximum allowed fraction of DM in the Milky 
Way contained in UCMHs, as a function of k and the UCMH 
mass MucMH- Here we show limits derived in this paper from 
Fermi-LAT searches for individual and diffuse DM sources. 
The UCMH mass is related to the mass contained inside the 
horizon when mode k enters by Eq. Q . All limits correspond 
to a 95% CL. Limits from searches for individual minihalos 
are based on non-observation of point or extended DM sources 
during one year of operation in all-sky survey mode. 



SUSY [I^H]^ as described in [H]. 

We assume 100% annihilation of WIMPs into bb pairs, 
a WIMP mass of = 1 TeV and an effective annihi- 
lation cross-section of {av) = 3 x 10~^^ cm^^ s~^. These 
are fairly conservative choices as far as gamma-ray yields 
go. Heavier WIMPs give lower fluxes, and WIMP masses 
considerably higher than 1 TeV are extremely challenging 
to obtain if one hopes for an associated natural solution 
to the gauge hierarchy problem. Our annihilation cross- 
section is the canonical, unboosted s-wave value implied 
by the relic density of DM under the assumption of ther- 
mal production. The bb final state gives rise to a rela- 
tively soft continuum spectrum dominated by pion decay; 
significant yields into final states with lower gamma-ray 
yields, such as /^^/i^, typically only arise in WIMP mod- 
els engineered to explain cosmic- ray excesses, and are in 
those cases accompanied by a corresponding boost factor 
in the annihilation rate. Taking both these effects into 
account, integrated UCMH gamma-ray fluxes from mod- 
els annihilating into /i^/i^ are not enormously different 
to those arising from unboosted annihilation into bb [IS] • 

Reading from right to left in Fig. [T] the Galac- 



tic gamma-ray source limit strengthens with increasing 
UCMH mass as UCMHs become brighter and more of 
the Galaxy is contained within dobs- At a mass of 
^7 X 10'^ M0, the limit is strongest, and UCMHs must 
constitute less than about 4x 10^'' of all DM in the Milky 
Way. 

At masses above ~10^ Mq, the value of /max given by 
Eq. (24) corresponds to less than three UCMHs in the 
entire Milky Way. At this point, Eq. (24) breaks down 



due to low-number statistics; the Milky Way may sim- 
ply contain zero UCMHs of a given mass in this case 
purely by chance, even though they are cosmologically 
more abundant than Eq. (24) would suggest. For larger 



masses, we obtain constraints by assuming that on aver- 
age there are at most 3 UCMHs of the mass in question 
per Milky-way sized halo. For zero observed UCMHs 
of a given mass, the CL with which we exclude a model 
that gives on average n such UCMHs per Milky- way sized 
halo is x[l — exp(— n)], with x the confidence level of the 
observation as in Eqs. ([24 26); choosing n = 3 makes 
this is a 95% confidence exclusion. From the minimum 
at intermediate UCMH mass, larger masses lead to less 
stringent limits on /, as three UCMHs progressively oc- 
cupy a larger fraction of the mass of the Milky Way halo. 

The limiting behaviour in Fig. [T] at very large and very 
small masses is simple to understand. For large UCMH 
masses, M^^^^ approaches the Milky Way mass, and 
because the limit here is given by the assumption that 
no more than three UCMHs of a given mass exist in 
each Milky- Way size halo, /max approaches one at ex- 



actly M?, 



UCMH 



For small UCMH masses, M, 



eventually shrinks 



^ UCMH routines will be included in a future public release of 
DarkSUSY. 



to such an extent that Eqs. (24) and (26) become greater 
than one. Common sense of course dictates that UCMHs 
cannot make up more than 100% of the Milky Way mass, 
but in our formalism this knowledge does not place any 
limit on the size of cosmological perturbations at such 
large k. This is because the Fermi limit is 'saturated' 
with respect to perturbations of this size; the pertur- 
bations can be arbitrarily large and still (by definition) 
unable to cause over 100% of the mass of the Milky Way 
to reside in UCMHs. For this reason, we do not give any 
limits for wavenumbers where /max — 1- For searches 
for individual Galactic gamma-ray sources, this corre- 
sponds to UCMH masses below ^10^^ Mq. In principle, 
one could obtain some bounds for larger wavenumbers 
by relaxing Eq. ([2]), and deriving direct limits on /3 as a 
function of M^, in cases where large amplitude density 
fluctuations would result in / = 1 before z = 10. This 
would be rather brave however, as it is not known to what 
extent the radial infall and absolute survival approxima- 
tions, which we rely on here, should be violated close to 
/=1. 

Such masses are already well into the regime in which 
kinetic coupling of DM might be expected to wash out 
structures such as UCMHs anyway. Indeed, this is 
an important general caveat at low masses: depending 
upon the specific particle DM candidate, the resultant 
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mass cutoff for UCMHs can be many orders of magni- 
tude larger than the smallest masses we consider here. 
We urge the reader to remember that the limits we 
present do not apply below the WIMP kinetic decou- 
pling threshold, and that this threshold should be calcu- 
lated on a per-model basis S74j. For e.g. neutralino DM, 
the corresponding size of the smallest DM halos at the 
time of equality may in principle be anything between 
M, = 10"" Mo and Mi = lO'^ M© [H], corresponding 
to fc^ax ^ — 10^ Mpc~"^ (although it is often assumed 
to be -10-*^ A/q, as for a 100 GeV pure Bino [TUgHTTT] !. 



C. Limits from gamma rays: individual 
extragalactic sources 

It is straightforward to extend the Galactic source 
analysis to sources residing outside the Milky Way. The 
approximate expression for /max close to the Earth, 
Eq. (25), in fact holds in this case, as the approxima- 



tion relies only on the probed volume being much smaller 
than that in which the number of UCMHs was assumed 
to be fixed. ^ In this case, the mass contained within dobs 
is given by 



MdKd^t,, = — (dobs 



drnax,Mw) 



Pc^x + 



MW ■ 



(28) 



where p^^x local cosmological background DM 

density. We show the resulting values of /max in Fig. [T] as 
a dashed purple line. As expected, the extragalactic lim- 
its provide increased sensitivity in the large-mass region, 
where the Galactic limits are set by the n — 3 condition. 
Having constructed both the Galactic and extragalactic 
limits as 95% CL upper limits, we see that the extragalac- 
tic curve tracks the Galactic one where dobs is comparable 
to -Rmw- At smaller fc, the extragalactic limit turns over 
and becomes stronger again, as unlike the Galactic DM 
profile, the cosmological background density is approx- 
imately constant with increasing d, providing a volume 
boost to the sensitivity for increasing UCMH mass. 



likely to be contaminated by astrophysical sources like 
pulsars and supernova remnants. The flux one would ex- 
pect from unresolved Galactic UCMHs in this direction, 
integrated in energy but differential in solid angle fi, is 



dJ", 



diff 



/ 



UCMH "'0 



Px{d)T{d)cfdd, (29) 



with d„ 



MW 



Here, i?MW is the virial radius 



of the Milky Way and tq the distance of the Sun from 
the Galactic centre; we take i?MW = 305 kpc |106| and 



ro — 8 kpc. Eq. (29) holds whether one refers to the raw 



flux F{d) or the energy flux F^{d), which differs from 
the raw flux only by an additional weighting factor of E 
inside the energy integral of Eq. (27). 



To obtain conservative upper limits on /, we compare 
this prediction to the total diffuse gamma-ray energy flux 
observed by Fermi at the Galactic poles. Referring to e.g. 
Fig. 1 in the 2-year LAT bright source catalogue |112| . 



this is about 1 x 10"''' GeVcm"^ s^^ sr~^ We derive 95% 
CL upper limits on / by demanding that the diffuse en- 
ergy flux from UCMHs in this direction does not exceed 
the total measured value times 1.2, which corresponds to 
an upwards variation of two times (Jeff = 10%, the maxi- 
mum systematic error in the LAT effective area over its 
energy range |112j . We implicitly assume in this proce- 
dure that statistical variation of the polar flux over the 2- 
year observing period is below the level of the systematic 
error coming from the LAT effective area, a reasonable 
approximation for our purposes. 

From the upper limits on / presented in Fig. [TJ we are 
able to determine upper limits on the mass variance of 
perturbations using Eqs. (15) and (21), assuming that 



the fractional UCMH content of the Milky Way is in- 
deed indicative of DM in the rest of the Universe. We 
use Brent's Method^ to numerically determine the value 
cr^ jj = cTmax corresponding to /max, and show the resul- 
tant upper limits on cr^ jj in Fig. |2| For this figure, we 



have combined the three different Fermi-LAT gamma- 
ray limits just discussed into a single best limit. 



D. Limits from gamma rays: Galactic diffuse 
emission 

For very small UCMH masses, where dobs is small but 
UCMHs are potentially very numerous, one might expect 
the diffuse gamma-ray background to provide a stronger 
limit on the UCMH fraction than searches for individual 
sources. In this case, a simple and robust limit can be 
obtained by considering the gamma-ray flux from a direc- 
tion perpendicular to the Galactic disk, as this is least 



E. Limits from reionisation 

We also plot in Fig. the upper limit on cr^ jj based on 
the limit on / at the ume of matter-radiation equality, 
calculated by Zhang [73]: 

/cq < 10-2(TOx/100GeV). (30) 

This limit comes from the impact of DM annihilation 
in UCMHs upon reionisation and the integrated opti- 
cal depth of the CMB at z < 30 (as seen by WMAP5 



Formally, one can replace Mmw in Eq- | |24| with e.g. the mass 

of the local group, cluster or supercluster, and make the same ^ Refer to standard numerical texts such as [113| for a description 
Taylor expansion as required to obtain Eq. \2b\. of all algorithms mentioned in this paper. 
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FIG. 2: Upper limits on the mass variance (J^^u horizon 
entry {aH = k), implied by the present-day UCMH abun- 
dance limits presented in Fig.[l] as well as the limits on / at 
the time of matter-radiation equality derived from the CMB, 
see Eq. (301. This later limit refers to the impact of UCMHs 
upon reionisation [73]; larger values of cr^ h correspond to a 
UCMH fraction that speeds up reionisation to the point where 
the integrated optical depth of the CMB (tg) is not consistent 
with the value measured by WMAP5 [114] . 



tion contained in massive compact halo objects (MA- 
CHOs) have been provided by EROS and OGLE 

(see e.g. [HI HH), on scales of A: - lO'' Mpc^i to 
k - 10^ Mpc-i. 

UCMHs are in fact just a non-baryonic variant of MA- 
CHOs. Naively applying these limits, which are very 
roughly of the order of /max ~ 0-lj would provide con- 
straints on the primordial spectrum of density fluctua- 
tions (see Section IV ) that are only slightly weaker than 



the constraints we derive from the gamma-ray limits. 
More importantly, these limits would not rely on the 
WIMP hypothesis and thus be completely independent 
of the nature of DM. 

As it turns out, however, UCMHs cannot be expected 
to have shown up in the EROS or OGLE data analy- 
sis because they are not sufficiently point-like. Indeed, 
their mass is much larger than the mass contained inside 
their Einstein radius, which significantly reduces the ex- 
pected magnification in lensing events. For this reason, 
we do not include any limits from gravitational lensing 
in our analysis. We note however that the light curve 
expected from UCMH lensing could help to detect such 
events even against a large background of non-MACHO 
events |48j . an effect which may produce interesting lim- 
its for future photometric lensing searches with, e.g., Ke- 
pler [125J . Another very promising way to constrain the 
abundance of UCMHs is to use astrometric microlens- 
ing, which involves searching for small changes in the 
apparent position of background stars with future planet- 
hunting facilities (126i il27] . 



[114) ). For larger /oq, reionisation takes place at higher 
z, introducing free electrons to the intergalactic medium 
at earlier times and increasing the Thomson-scattering 
optical depth to the surface of last scattering (r,,). We 
implemented this limit with the same DM mass of 1 TeV 
discussed above. According to Eqs. Q, ([8| and (27), 



the number of ionising photons produced by a sing. 
UCMH scales linearly with its mass, as 



oc J-(rf)d2 



r^p^{r)dr cx M°cmh- (31) 



For a fixed / the number of UCMHs scales inversely with 
their mass, so we see that the limit in Eq. ( 30 ) should be 
independent of the UCMH mass. This limit is weaker 
than those from gamma rays. We do not give this limit 
in Fig. [l] as such a UCMH fraction at equality would 
grow to become simply f — 1 today (i.e. the bound 
saturates before z = 0). 



IV. CONSTRAINTS ON THE PRIMORDIAL 
POWER SPECTRUM 

A. Constraints from UCMHs 

Using the limits on the mass variance of Gaussian per- 
turbations presented in Fig. [2j we can now begin to con- 
strain the primordial power spectrum. This requires an 
explicit relationship between cr^ ^ and the power spec- 
trum itself. We consider three broad models for the spec- 
trum of perturbations: 

1. A scale-free spectrum 'P{k) (x fc"~^ parameterised 
by a spectral index n. 

2. A scale-free spectrum with a step, resulting in in- 
creased power on small scales. This spectrum is 
parameterised by the spectral index n, as well as 
the height p and location kg of the step: 



F. Limits from gravitational lensing 



r{k) oc X 



for k < ks 
for k > k^ 



(32) 



Various gravitational microlensing experiments look- 
ing towards the Large and Small Magellanic clouds have 
been carried out in recent times (115H121] . The strongest 
resulting limits to date on the Milky Way halo mass frac- 



3. A non-parametric generalised spectrum, assumed 
to be locally scale- free in fc-space, i.e. Vik) oc k'^~^, 
but in general allowing for a different normalisation 
for very different values of k. 
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FIG. 3: Upper limit on the spectral index of the primordial 
power spectrum from gamma-ray searches for UCMHs, and 
the impact of UCMHs on reionisation. These limits assume 
5^ oc k"~^ , and take into account only the constraints on cr^,H 
given in Fig.[2j for wave numbers smaller than k. For compar- 
ison, we also show the resulting gamma-ray constraint if we 
were to assume 5™'" — 10"'^ (improved upon in Appendix A| 
and use the over-simplified calculation of a"x,H US] (corrected 
in Appendix |B|) . 



In the first two cases, we normalize our spectra to the 
WMAP data in the same way as done in Ref. |128j (which 
is fully consistent with, and actually only marginally im- 
proves, the normalization measured by COBE '50'). For 
a more detailed description, and a derivation of the rela- 
tionships between cr^ jj and the parameters (or running 
amplitude, in the generalised case) of each spectrum, see 
Appendix [B] 

Using these expressions, we translate the upper limits 
on cr^ H in 
each spectr 

index n of the scale-free spectrum, derived from UCMHs 
of different masses, corresponding to different values of 
k. Using such a spectrum, the most relevant limit is es- 
sentially the lowest one obtained at any scale. Using this 
criterion, our limits on the UCMH abundance constrain 
the spectral index to be n < 1.17. 

In the case where the spectrum contains a step, we de- 
rive upper limits on the step size p as a function of its lo- 
cation fcg (Fig. |4j left panel), assuming the spectral index 
n = 0.968 observed on large scales by WMAP 7 [51]. We 
also show how these limits change if n is allowed to vary 
within the la range of the WMAP7 measurements, plot- 
ting the resultant regions as shaded bands. Depending 
upon the true value of n, the non-observation of UCMHs 
limits the size of a step in the power spectrum at any 
scale larger than fc = 3 x lO^Mpc"^ to be less than a 



Fig. 2 into upper limits on the parameters of 
rum. Tn Fig. [3j we give limits on the spectral 



factor of 10-12, quite a severe constraint. 

For each value of kg, the limit on p is based on a par- 
ticular scale k (and therefore one specific M^^^-^). For 
each ks, we chose the optimal value of k by minimising 
the resultant Pmax over all k using Brent's Algorithm. 
As is to be expected, optimal limits always come from 
length scales smaller than the scale of the break, as this 
is the region affected by the step. Because the limits on 
CT^ H from gamma rays and reionisation are (albeit only 
approximately, in the case of gamma rays) monotonically 
decreasing functions of fc, the optimal k for placing limits 
on p is always the smallest scale considered valid. In our 
case, the dominating scale is A; = 3 x lO^Mpc"^, which 
corresponds to our minimal assumed free-streaming mass 
cutoff, Mi = 5x 10-12 M0 (M^cMH ~ 10"^ Mq). The 
limits in the left panel of Fig.[4]are therefore entirely dom- 
inated by the free-streaming scale, and should be treated 
with caution, as they will weaken for DM models that 
lead to larger minimal halo masses. In the right panel of 
Fig. |4j we give an example of how the limits on p for a 
spectrum with a step at kg = lO^Mpc"'^ weaken if the 
kinetic decoupling scale is varied. Here we consider the 
range Af,; = 3 x 10^"' — 5 x lO^^^ Mq, corresponding to 
allowed values within an indicative set of minimal super- 
symmetric standard model (MSSM) benchmark points 
[74] . As expected, larger minihalo mass cutoffs weaken 
the corresponding limits on the size of a step in the pri- 
mordial spectrum; the exception to this behaviour occurs 
for gamma-ray limits below ~3 x lO^Mpc"^, where the 
limits on cr^ jj in Fig. j2jare not monotonically decreasing 
with k. 

In the left panel of Fig. [5] we give limits upon the am- 
plitude of curvature perturbations V-jz at small scales, 
using a non-parametric description of the spectrum. The 
corresponding limits on the amplitude of physical den- 
sity perturbations Vs can be obtained by simply mul- 
tiplying the curvature perturbation limit by a factor of 
0.191 (Eq. B23 ). For comparison, we also show the rather 



different limits on curvature perturbations presented in 
Ref. [55], derived from Fermi searches for DM annihila- 
tion. We attribute the differences in the two limits to the 
scale-dependent calculation of Smin we perform here, our 
improved statistical treatment of the limits coming from 
Fermi, the inclusion of diffuse gamma-ray limits, and our 
corrected calculation of the mass variance of perturba- 
tions relative to Ref. [55] , 

We also show in the right panel of Fig. [5] the depen- 
dence of our limits upon the WIMP mass, and the latest 
allowed redshift of UCMH collapse. As argued earlier, 
our choices clearly are indeed very conservative. Smaller 
WIMP masses and, to a much larger extent, later red- 
shifts of collapse result in significantly strengthened lim- 
its - simply because in the latter case much smaller initial 
density perturbations would have time to collapse, see 
Appendix [A] and in particular Fig. [7j The minimum red- 
shift of collapse could actually quite defensively be made 
somewhat smaller, without invalidating our earlier argu- 
ments about the UCMH survival probability and viola- 
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FIG. 4: Left: Constraint on the allowed height p of a step in the primordial power spectrum from gamma-ray searches for 
UCMHs and impacts on reionisation, as a function of the location ks of such a step. Here p refers to the dimensionless ratio 
of the power at the wavenumbers immediately above and below ks- Our central curves assume the spectral index n = 0.968 
obtained from WMAP7 observations of large scales [ST], and shaded regions correspond to the 68% CL for this measurement 
(An = 0.012). Right: Variation of the gamma-ray and reionisation constraints on p with the kinetic decoupling scale of DM. 
These two limits in particular are sensitive to the cutoff in the DM halo mass function, as the strongest limits (as shown in the 
left panel) come from the smallest viable UCMHs, for all ks- 



tion of the radial infall approximation. In the interests of 
producing as robust limits as possible, we use Zc = 1000 
as our canonical value, but it is worth noting that the 
substantial improvement in limits with smaller Zc bodes 
well for the potential for future detection of UCMHs. 



B. Comparison with existing constraints 

On large scales, the primordial power spectrum of den- 
sity fluctuations has been measured with high precision 
mainly by CMB experiments (as the most powerful sin- 
gle source of information about the primordial fluctua- 
tions), large scale galaxy surveys (large scale structure; 
LSS) and weak gravitational lensing observations. How- 
ever, the scales probed by such cosmological measure- 
ments constitute only a relatively small part of the entire 
spectrum, namely scales between k ~ 10""* Mpc~^ and 
A: ~ 1 Mpc~^. This range has been extended to smaller 
scales (fc 3 Mpc~^) by other astrophysical measure- 
ments that probe later epochs in the evolution of the 
Universe, such as the Lyman-a forest. For the rest of the 
spectrum, i.e. on scales smaller than ~ 3 Mpc~^, con- 
straints have been provided mainly from non-observation 
of PBHs. Although PBHs constrain the power spectrum 
over a very wide range of small scales (from k ^ 10^^ 
Mpc~^ to fc ~ lO^'^ Mpc~^), the constraints are different 
from those provided by the aforementioned cosmologi- 
cal probes in that 1) they are upper limits rather than 



positive measurements and 2) are much weaker (the up- 
per limits are many orders of magnitude larger than the 
cosmological constraints on large scales). 

A number of different techniques have been used to 
constrain the primordial power spectrum and its proper- 
ties with cosmological observations. These include top- 
down approaches, where specific theoretical models of 
the primordial fluctuations or the inflaton potential are 
fit to the data, and bottom-up attempts to reconstruct 
the shape of the spectrum from data with no such as- 
sumptions. Such reconstruction techniques consist of 
simple binning techniques I 129H138] . principal compo- 
nent analysis 139 , methods of direct inversion (1401 - 
I154j and minimally-parametric reconstructions based on 
cross-validation technique [1551 1156j . Results of these re- 
construction procedures, although consistent with each 
other, are in general not identical even in cases where 
the same observational datasets are employed. 

One of our objectives in this paper is to present a 
comprehensive compilation of the best constraints on 
the power spectrum at all scales, including those from 
UCMHs (as discussed in Sec. IV A). We therefore select 



the strongest available constraints provided by the latest 
analyses of different cosmological data; these come from 
Refs. dSl [m [ISS], and draw primarily upon CMB, 
LSS and Lyman-a data, though some small additional 
constraining power is derived from measurements of pri- 
mordial nucleosynthesis, supernovae and the Hubble con- 
stant (refer to Ref. }153j for details). We combine these 
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FIG. 5: Left: 95% CL upper limits on the amplitude of primordial curvature perturbations Vn (for a non-parametric, generalised 
spectrum) allowed by gamma-ray searches for UCMHs and impacts on reionisation. For comparison, we show the previous 
limits from Fermi non-observation of UCMHs derived in Ref. [55], based on a simplified treatment of the statistics of non- 
detection, mass variance and minimum density contrast required to form a UCMH. Corresponding constraints on the generalised 
amplitude of primordial density perturbations can be obtained by multiplying these limits by a factor of 0.191, see Eq. (B23l. 



Right: The variation of the gamma-ray limit on Vn with WIMP mass and the redshift of UCMH collapse, showing the impact 
of less conservative (but entirely plausible) choices for these parameters than our canonical = 1 TeV, Zc — 1000. 



constraints into a single Icr band (comprised of the best 
available upper and lower limits at each k). We plot the 
resultant constraint band in Fig. |6j along with our own 
results on small scales from reionisation and gamma-ray 
searches for individual UCMHs and Galactic diffuse DM 
sources. 

We also show the current strongest upper limits on the 
power spectrum derived from PBHs [33], based on their 
present-day gravitational influence (for k < lO^^Mpc"^) 
and the products of their conjectured [157] evaporation 
(for k > 10^^ Mpc~^). We do not show limits above k > 
10^^ Mpc^^, as such constraints rely on model-dependent 
assumptions about new (quasi-)stable elementary par- 
ticles that often appear in extensions of the standard 
model. At even smaller scales, Ip^ k > lO^^Mpc""'^, 
the situation is even more uncertain: one must assume 
that the evaporation leaves a hypothetical Planck-size 
relic in order to place any further limits. 

We do not make any attempt to harmonise the CLs 
with which we state limits from different sources, as we 
do not have access to the full likelihood functions for 
any of the reported results. Let us also repeat a general 
word of caution for these kinds of limits: even though we 
present them here in a model-independent way, such lim- 
its always depend, to a certain extent, on the assumed 
spectrum of the density perturbations. This should be 
kept in mind when comparing predictions from, e.g., in- 
flationary models to what is shown in Fig. [6] (see Ap- 
pendix |B] for how to treat spectra that deviate from the 



locally scale-free spectrum that we assumed here). 

V. CONCLUSIONS 

In ACDM cosmology, the first gravitationally bound 
objects typically form at redshifts considerably smaller 
than z ^ 100. Very rare fluctuations, on the other hand, 
would collapse as early as the first stages of matter domi- 
nation, forming UCMHs. By the time standard structure 
formation starts, these objects would have already devel- 
oped a dense and highly concentrated core, which would 
be essentially unaffected by tidal interactions during sub- 
sequent cosmological evolution. 

Arguably, the most compelling class of candidates for 
the nature of DM are WIMPs. Because UCMHs have 
survival probabilities close to unity, the non-observation 
of individual gamma-ray sources by Fermi, as well as 
measurements of the diffuse gamma-ray background, can 
then be used to place rather strong limits on the power 
spectrum of primordial density perturbations. Potential 
changes in the reionisation history of the Universe, which 
would leave a visible imprint in the cosmic microwave 
background, can also provide some constraints. 

We have provided and discussed in detail all necessary 
formulae to calculate these constraints for any functional 
form of the primordial spectrum of density fluctuations. 
As a possible application, we have constrained the spec- 
tral index of an assumed featureless power-law spectrum 
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FIG. 6: Constraints on the allowed amplitude of primordial density (curvature) perturbations Vs {Viz) at all scales. Here we 
give the combined best measurements of the power spectrum on large scales from the CMB, large scale structure, Lyman-a 
observations and other cosmological probes 152, 1531 1156] . We also plot upper limits from gamma-ray and reionisation/CMB 
searches for UCMHs, and primordial black holes |45] . For ease of reference, we also show the range of possible DM kinetic 
decoupling scales for some indicative WIMPs [74]; for a particle model with a kinetic decoupling scale Aikd, limits do not apply 
at fc > fcxD. Note that for modes entering the horizon during matter domination, Vs (but not Vn) should be multiplied by a 
further factor of 0.81. 



to be n < 1.17. Since large-scale observations actually 
put much stronger limits on the spectral index, we have 
also considered the case of n = 0.968 ± 0.012, as ob- 
tained by WMAP observations, and constrained the al- 
lowed additional power below some small scale fcg to be at 
most a factor of ~ 10-12 (assuming a step-like enhance- 
ment in the spectrum). As a third example, we have 
obtained quasi-model-independent limits, of the order of 
Vn < 10~^, on perturbation spectra that can at least 
locally be well described by a power law. We would like 
to stress, however, that it is intrinsically impossible to 
constrain primordial density fluctuations in a completely 
model-independent way; one thus has to re-derive such 
limits for any particular model of, e.g., inflation which 
produces a spectrum that does not fall into one of these 
classes. Here, we have provided all the necessary tools to 
do so. 

We have mentioned that present gravitational lens- 
ing data cannot be used to constrain the abundance of 
UCMHs - essentially because they are simply not point- 
like enough, even in view of their highly dense and con- 
centrated cores. Future missions making use of the light- 
curve shape in lensing events, however, are likely to probe 
or constrain their existence. This would be quite remark- 
able as it would allow us to put limits on the power spec- 
trum without relying on the WIMP hypothesis for DM. 
Most of our formalism is readily extended, or can in fact 
be directly applied to, such constraints arising from grav- 
itational microlensing. 

Finally, we have compiled an extensive list of the most 



stringent limits on V-jiik) that currently exist in the lit- 
erature for the whole range of accessible scales, from the 
horizon size today down to scales some 23 orders of mag- 
nitude smaller. Direct and indirect observations of the 
matter distribution on large scales - in particular galaxy 
surveys and CMB observations - constrain the power 
spectrum to be Vn{k) ~ 10~^ on scales larger than about 
1 Mpc. On sub-Mpc scales, on the other hand, only upper 
limits exist. From the non-observation of PBH-related 
effects, one can infer Vtz < 10"^ — 10"^ on all scales 
that we consider here. UCMHs are much more abun- 
dant and thus result in considerably stronger constraints, 
T^n ^ 10"^, down to the smallest scale at which DM is 
expected to cluster (this depends on the nature of the 
DM; for typical WIMPs like neutralino DM, e.g., it falls 
into the range fc^^^ - 8 x 10"* - 3 x lO^Mpc"^). 



It is worth recalling that the observational evidence 
for a simple, nearly Harrison-Zel'dovich spectrum of den- 
sity fluctuations is obtained by probing a relatively small 
range of rather large scales. The limits we have provided 
here will thus be very useful in constraining any model of 
e.g. inflation, or phase transitions in the early Universe, 
that predicts deviations from the most simple case and 
which would result in more power on small scales. 
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Appendix A: Minimal density contrast for UCMH 
formation 



In their original paper, Ricotti & Gould [JS] used an 
order-of-magnitude estimate to argue that a region of co- 
moving size R should have an average over-density of (5 > 
jmin _ ]^Q-3^ g^i- ^Yie time of horizon crossing,^ in order 
to gravitationally collapse to a UCMH before a redshift 
of 2; ~ 1000. This is also the value that has been used 
in subsequent work, see e.g. Refs. [IHl [SSJ [73] . However, 
since (5™™ enters exponentially in the expression for the 
abundance of UCMHs, it is important to have a more 
accurate estimate of this quantity. Note, in particular, 
that one should expect it to be scale-dependent rather 
than just a constant since DM perturbations continue to 
grow after they enter into the horizon (even though this 
growth is only logarithmic for times long before matter- 
radiation equality |158] ). 

To derive the minimal density contrast required for 
UCMH formation, we will in the following rely on a 
simplified description for the collapse of over-dense re- 
gions that was originally introduced in |159| and is now 
widely being used. The basic approximation is to re- 
strict oneself to the case of spherical regions of uniform 
density p{t) = p{t) + Sp{t) that are embedded in a back- 
ground with density p{t). During matter domination. 



In Fourier space, the time tk of horizon crossing is convention- 



ally defined as a,H\^ 



k, where k is the wavenumber of the 



perturbation. In real space, we define it as aR = H~^\ this 
amounts to saying that we associate a wavenumber kft = R^^ 
to a perturbation of co-moving size R. 



the full (non-linear) evolution of the over-density can in 
that case be solved exactly; p{t) starts off by decreasing 
more slowly than the background density p{t), reaches a 
minimum after a finite time tc/'i and then increases until 
it becomes infinitely large at t = tc- In this approach, 
the time tc is thus identified as the time where the re- 
gion has fully collapsed - and which at the same time 
indicates the breakdown of this simplified description (in 
a more realistic description, dynamic relaxation and an- 
gular momentum conservation would of course prevent 
the collapse to a point-like, singular object). In the lin- 
ear theory, on the other hand, the over-density would by 
that time have grown to 



5p{t 



cjlin 



P{tc 



2/3 



1.686. 



(Al) 



This relation thus allows us to use the linear theory of 
perturbations to calculate the time of collapse (at least 
during matter domination) - which is quite remarkable 
since the perturbations of course enter the non-linear 
regime already quite some time before the actual col- 
lapse. 

Unfortunately, even in linear theory the system of 
equations governing the evolution of density contrasts 
around the time of matter-radiation equality is rather 
complicated and cannot be solved analytically. However, 
a very accurate fit to a numerical solution for the density 
contrast in DM fluctuations for t > t^q can be found, e.g., 
in Ref. [7j: 

S^{k,t)^^r{n)n'{k), (A2) 

where k is the co-moving wave number of the perturba- 
tion and K its rescaled, dimension-less version: 



/2k 



The fitting function T{k) is given by 
ln[l + (0.124k)2] 



r(^) 



(0.124k)2 

1 + (1.257k)2 + (0.4452k)4 
1 + (1.606k;)2 -f- (0.8568k)4 



(0.2197k)^ 
(0.3927K)f 



(A3) 



(A4) 



1/2 



The above expression for d^. is given in synchronous 
gauge, but since we are interested in scales that are 
much smaller than the horizon at the time of collapse, 
the choice of gauge does not actually matter. The nor- 
malisation is chosen such that, for adiabatic fluctuations, 
TZ'^ gives the value of the initial curvature perturbation 
(which is time-independent on scales much larger than 
the horizon). 

By equating Eqs. (Al) and (A2), we can thus derive 



the minimal value of 7?," that is required so that the per- 
turbation in the DM component collapses before a given 
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redshift Zn < Zr,, 



a 



(A5) 



While TZ'^ is the actual fundamental physical quantity 
describing the strength of adiabatic fluctuations - pre- 
dicted, e.g., in theories of inflation - one may instead 
also consider the somewhat more intuitive value of the 
DM density contrast at the time a fluctuation enters into 
the horizon - which in our case would be during the 
radiation-dominated era. In contrast to 7?.°, however, the 
density contrast is a gauge-dependent quantity (which 
becomes numerically quite relevant for scales k < aH). 
In the following, we will choose the co-moving ( or "total 
matter") gauge, where the rest frame is that of the total 
energy density fluctuations; this is the gauge that corre- 
sponds to the initial conditions adopted in the treatment 
of the collapse outlined above (and also is typically used, 
e.g., for the calculation of PBH formation). 

During radiation domination, the evolution of all per- 
turbed quantities like 6-^ can be solved analytically. To 
convert the results in synchronous gauge given in Ref. [7; 
to the total matter frame, one has to perform a gauge 
transformation under which 



(A6) 



where p is the time derivative of the mean density and 
is the scalar part of the velocity components in the 
stress-energy tensor in synchronous gauge. ^ The result 
can be written as 



6,{k,t) = N,0^T,{e)n''(k), 



where 



2kt 



1 k 



(A7) 



(A8) 



and I = X)*^ for the DM and radiation component, re- 
spectively. For adiabatic fluctuations we have Nr = 4/3 



and 



1. The transfer functions introduced here 



satisfy Ti(0) = 1 and are given by 



6 



ln(6l) + 7E 



Ci(^) + xJo(^) 



(A9) 
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FIG. 7: Minimal density contrast in the DM component, at 
the time of horizon entry (fc = aH), required to form a UCMH 
before redshift Zc = 1000, 500, or 200. 



where 7e is the Euler-Mascheroni constant, Ci the cosine 
integral function and jp i a re spherical Bessel functions 
of the first kind. Eq. (AlO) is somewhat painful to im- 



plement numerically in certain good solution is 

to modify an existing code for computing Ci(a;) (as found 
in e.g. iUSj) to instead return x~^[Ci{x) — In a; — 7e]. 

We have now collected all the pieces needed to express 
the minimal density contrast in the DM component, at 
the time of horizon entry tk of a scale k, for a perturba- 
tion to collapse and form a UC MH before a redshift Zc- 
Combining Eqs. (A5) and (A7), we have 



{k,tk) 



a 



2/3 



T^je = 1/73) 

k^TiK) 



(All) 

This function is plotted in Fig. [7] for a few selected values 
of Zc- Note that after decoupling at Zdcc ^ 1000, baryons 
would also start to gravitationally collapse and thus sig- 
nificantly contribute to the over-density - which has not 
been taken into account here. In that sense, our estimate 
for 5™'" is rather conservative for redshifts z <§: Zdec- 



This transformation behaviour follows from the fact that both 
synchronous and total matter gauge have vanishing metric com- 
ponents goi and that in co-moving gauge we have Su = 0. Note 
that the latter condition follows from 5u being defined as part 
of the stress-energy tensor (as in Ref. |7|) — which is a somewhat 
different definition compared to other examples from the litera- 
ture I160| : as a result, the change of velocity perturbations under 
gauge transformations does not take the same form as in these 
references, either. Of course, the final results are unaffected. 



Appendix B: Correct normalisation of power 
spectrum and mass variance 

In this Appendix, we review in detail how to express 
the mass variance, i.e. the r.m.s. over-density in a given 
region of space, in terms of the super-horizon spectrum 
of density or curvature fluctuations provided by inflation. 

Assuming Gaussian statistics for the primordial den- 
sity fluctuations, the probability (density) to find an av- 
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erage density contrast 5 in a spherical region of size R is 
given by 



PRiS) 



exp 



2a2(i?) 



(Bl) 



where the mass variance a{R) is computed by convolving 
the power spectrum with a top-hat window function: 



a^R) = f 
Jo 



w^^{kR)rs{k) 



dk 



In the above expression, Vs{k) is defined by 



n(fc)'5(k-k') 



(B2) 



(B3) 



and W^ii(x) — 3ji{x)/x = 3x~^ (sin a: — a; cos a;) denotes 
the Fourier transform of the 3D top-hat window function, 
with X = kR. 

The above description is somewhat complicated by the 
fact that density perturbations evolve with the expansion 
of the Universe, which means that also the power spec- 
trum is time-dependent. In total matter (as well as in 
synchronous) gauge, the quantity^" 



5Uk,t)^^^Vs{k). 



(B4) 



however, is time-independent on super-horizon scales 
(k <^ aH), with a numerical value that is very close to 
the value at the time when mode k crosses the horizon. 
It is thus illustrative to separate the power spectrum into 
a part that describes the primordial fluctuation spectrum 
on super-horizon scales, as provided by inflation, and 
a part that encodes the evolution of the perturbations 
(mostly after they enter the horizon). This is usually 
done by introducing a transfer function 



T^{k,t) = 5l{k,t)/5l{k,U 



(B5) 



which satisfies T{k 0,t) — 1. In the above definition, 
ti denotes a time before the entrance of any scale k into 
the horizon, U < tk Vfc, so it should be taken to corre- 
spond to the time at the end of inflation. Note that the 
time-dependence of the transfer function only enters as 
the ratio of aH/k, so T can also be written as a func- 
tion of one variable only. For e.g. radiation domination, 
explicit expressions are given in Eqs. (A8 AlOl. The re- 
lation between S^i and the initial spectrum of (adiabatic) 
curvature perturbations during radiation domination is 



as can be explicitly verified with the help of Eqs. (A8 



A9). During matter domination, the right hand side of 
this expression has to be multiplied by (10/9)^ |161| . 

As explained in Appendix [S] the condition for UCMH 
formation can be stated in terms of the DM density con- 
trast at horizon crossing of a given scale. We should 
therefore evaluate the mass variance in Eq. ( |B2[ ) at the 
time when k^ = 1/ R = aH. With the shorthand nota- 
tion 



fe(ifc) = fe(fc,ifc) , 
this can be conveniently expressed as [101 03] 



4(i?) = ^2 {R)\ 



(B7) 



(B8) 



where 



a'{k)^5^^{tk) / xHl{xk,U)T\xk,tk)W^^{x)dx . 
Jo 

(B9) 

The important point to note here [101 HE] is that the re- 
lation between mass variance (Th and size of the pertur- 
bation at horizon crossing 6^ depends on the scale and, 
in principle, on the full cosmological evolution between 
the end of inflation and tk- 

Unlike for the case of PBHs, only the density contrast 
in the DM component will grow and eventually collapse 
to a UCMH; in the above expressions, we should thus 
place an index x where appropriate. Since S and only 
differ by a constant factor on super-horizon scales (at 
least for curvature perturbations), the main difference 
between a{k) and a^{k) is the transfer function that is 
used. This difference, however, is important because per- 
turbations in the DM component behave completely dif- 
ferently from those in other components once they enter 
the horizon, cf. Eqs. (A9-A10). The mass variance we 



are really interested in is therefore given by 



alikn)S^{tkJ 



(BIO) 



where S^{tkii) refers to the total energy fluctuations be- 
cause we later want to normalise it to the present-day 
density contrast observed in the CMB. For k ^ fceq = 
\/2Honjn/V^i 'v^e have 



9 
16 



dxx^W^T^{x) 



rn{xk)T^{e = x/V3) 
Vn{k) T2(0 = l 



3) 
(Bll) 



VTzik) 



S^iik, ti), 



(B6) 



The quantity 5h is directly related to what is also known as the 
peculiar gravitational potential, <1?^ = (97r^/2A;'^)(5g. 



1. Scale- free spectrum 

It is often assumed that the spectrum on super-horizon 
scales is of a scale-free form, implying also that 



SUtk) oc fc"- 



(B12) 
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during any epoch where the equation of state does not 
change. In this case, Eq. (Bll) simphfies to 



TH 



T2 

(x) ^ 



r2(i 



(B13) 



For a spectrum Uke in Eq. (B12) and a flat cosmol- 
ogy, the WMAP normalisation at the scale /cwmap = 
0.05Mpc~i is given by [HHI 

A„ = ia.7.in-^-P[(^- -)(-^-^4 + ^-»4-)], (B14) 

Vl + 0.53r 

where r is the tensor-to-scalar ratio. Assuming no grav- 
itational waves (i.e. r = 0) and taking into account that 
the normalisation of TZ leads to an effective suppression of 
5h by a factor of 9/10 during the transition to the matter 
dominated phase [161] . we thus have for ^ t^q- 



10\ 2 

^ j ^H(*fcwMAp) 

4.58 X io-i0e2.48(n-i) 



k 

Wmap 

k 



(B15) 



fcwMAP 



For non-zero values of r, in principle, this normalization 
changes slightly. From inflation in the slow-roll approx- 
imation, e.g., we expect r > (8/3) (1 — n) ~ 0.085 (as- 
suming the spectral index observed on large scales), cor- 
responding to a very modest decrease of by about 2%. 
Even the largest currently allowed value of r 0.3 [51] 
would decrease 5^ by only about 6%, leading to a cor- 
responding weakening of our limits on the height p of a 
step in the spectrum (see below) and barely affecting our 
limits on n. 

We now have all the ingredients to calculate cr^jj as 



in Eq. (BIO). Approximating the transition between ra- 
diation and matter domination to take place almost in- 
stantaneously, the last factor in the above expression can, 
furthermore, be expressed in terms of the total mass con- 
tained within the horizon as 





■ Mnito) ' 


(n-l)/3 




(n-l)/2 








_MH{tk)_ 





(B16) 

where the first term on the right-hand side can also be 
written as {1 + z^qY^^^ . 



2. Spectrum with a superimposed step 

The most simple phenomenological way to allow for 
more power below some small scale (k ^ fcoq) is to gen- 
eralise the scale-free case by introducing two parameters 
{ks,p) that describe the location and height of a step in 
the primordial spectrum (see e.g. [TOl 



"H.step 



1 for fc < kg 
for k > kg 



(B17) 



where S'^(tk) is given by Eq. (B15|. We can then again 
use Eq. (BIO) to calculate a^^-a by using 

9(p2 - 1) 



{k, ks,p, n) = a (n) + 



16 



(B18) 



1a "^^^"^,2(1/73) 



with a^{n) given by Eq. (B13 ). Note that a rather similar 



spectrum can arise if the infiaton potential has a jump 
in its first derivative [T2l. 



3. Generalised spectrum and curvature 
pert ur bat ions 

Instead of relating it to (5h, we can of course also ex- 
press the mass variance directly in terms of the curvature 
perturbation and write Eq. (BIO) as 



= IJ x^W^^ix)Vnix/R)T^ (x/^/3) dx. 

(B19) 

This integral is strongly dominated by contributions 
around x 1 (i.e. by wavenumbers k^ — 1/R). A gen- 
eral form for the spectrum of primordial curvature per- 
turbations can be obtained by discarding the assumption 
of a globally scale-free spectrum, and instead adopting 
the weaker assumption of local invariance only around 
the scale of interest k^ probed by (non-)observation of 
UCMHs of a certain mass (see e.g. [SS]). Given the dom- 
inance of /cr, in the integral above, this can be achieved 
by taking a power- law form for Vn{k), with the pivot 
point placed at /cr. 



k 



(B20) 



Here ^^(/cr) is the local slope of the power-law at fcpj. 
For e.g. slow-roll inflationary models, 0.9 < riR < 1.1 
is typical at scales measured by the CMB \16'S\ : for the 
considerably smaller scales we are mainly interested in 
here, however, it is much more difficult to make general 
statements about the expected value of n^j. For our con- 



straints we use nR = 1, for which evaluation of Eq. (B19 ) 
gives 

<h(^)/^k(%) = 0.907. (B21) 

For comparison, changing 71r to (2) would result in 

0. 388 (2.91) - which should serve as a warning that it 
is in general not possible to translate bounds on into 
bounds on V-ji in a completely model-independent way, 

1. e. without assuming anything about the form of the 
spectrum. 

Similarly, using Eqs. (B4 B6) we can express the 



equivalent total density power spectrum in terms of the 
initial curvature spectrum as 



Vs{k) 



2k y 



T^{k,t)Vn{k), 



(B22) 
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which at the time of horizon entry gives 



Vs{k)\t=t 



T^{9 = l/V3)Vn{k) = 0.191Vn{k). 

(B23) 
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